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AN ANALYSIS OF THE OPTIMIZATION OF A BEAM RIDER 
MISSILE SYSTEM 

By Marvin Shinbrot and Grace C. Carpenter 
SUMMARY 

A transfer function is derived for a beam rider missile guidance 
system -which is optimum when the target moves in a nonstationary way. 

The effects of acceleration limiting are considered and a discussion of 
the miss as a function of the various parameters which determine it is 
included. A form of design chart is presented which allows the immediate 
determination of the optimum under any set of target-missile conditions. 



INTRODUCTION 


With the recent progress in information theory and related techniques 
for optimization of systems operating in the presence of noise have come 
several applications of these optimization methods to the design of guided 
missile control systems (refs. 1 and 2). The results of these analyses 
have been transfer functions which specify the system which is optimum 
under the conditions assumed. The transfer functions thus obtained appear 
quite satisfactory, the mean-square miss distance associated with them 
being reasonably small and the functions themselves being not of an 
■unrealistic form. 

In order to arrive at their results, the authors of references 1 
and 2 used the classical Wiener theory (ref. 3)* As is well known. In 
order to apply the Wiener theory, it is necessary that the class of inputs 
to the desired system be stationary. Now, real targets may or may not 
maneuver in a nonstationary way. However this may be, one can easily find 
examples to which the Wiener theory as originally conceived does not apply: 
an example would be the case where the target maneuver consists of a step 
in acceleration. In such a case, if the Wiener theory alone were available 
as a tool to the designer, he probably would approximate the nonstationary 
maneuver by a stationary one. This implies that some improvement of such 
systems might be expected if more general target motions - ones involving 
no approximation - could be considered. 

A method which allows the optimum to be determined when the inputs 
are not stationary was presented in references h, 5, and 6. It is the 
purpose of this report to apply this method to the optimization of a beam 
rider control system. Our purpose here Is exemplary; although 311 optimum 


•* 



2 


NACA TN klk3 


transfer function is presented, we shall "be more interested here in 
indicating how the method can he applied to missile problems t han in 
specifying an optimum one. Thus, the point of view of the preceding para- 
graph will not he adhered to strictly. Although certain assumptions will 
he made which simplify the work, the situations described here will he more 
realistic than those of references 1 and 2. It should he stressed, how- 
ever, that the methods of references h, and 6 are sufficiently powerful 
that these assumptions can he eliminated. 

It might he of interest, while on the subject of assumptions, to 
discuss one which is not made. In all the previous works on missile opti- 
mization, it was assumed that the target and the missile move in one plane. 
One may he certain, however, that if the pilot of a target bomber knows 
this, he will do his best to assure that he does not remain in the same 
plane as the missile. Consequently, we shall not make this assumption: 
the plane in which it will he assumed the target moves will have no invar- 
iant relation to the missile’s initial position and the flight path of the 
missile will not he assumed to he planar at all. 

The paper begins with a discussion of the optimum system and the 
corresponding minimum error. It Is shown that this system can never 
actually achieve the indicated minimum miss, since in order to do so the 
level of acceleration required of the missile would he Impossibly high. 

The effect of introducing a side condition that the root-mean-square accel- 
eration shall not exceed a preassigned value Is then considered. By a new 
method another system is derived; this system has the properties that its 
rms acceleration remains below the given value while the miss corresponding 
to it is, for a fairly wide range of limiting accelerations, only slightly 
larger than the minimum miss. A brief discussion of this miss as a func- 
tion of the parameters which determine it is then given. 


SYMBOLS 


a vertical acceleration of target, ft/sec 2 

aj£ acceleration of missile, ft/sec 2 


E root-mean-square miss distance, ft 

HT N x +Ny+N z ' ' 

N x ,Ny,N z noise amplitudes (magnitudes of spectral densities at zero 

frequency) in the directions of the coordinate axes, ft 2 sec 


R 



Ir 


t 


time measured from instant of firing, sec 
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time at which target "begins maneuver, sec 


t* 


T 

V 


Xo^O^O 

x B^ y B ,Z B 


X M ,Z M 


"At, dimensionless time 

time at which missile velocity is so reduced, it can no longer 
capture target, sec 

horizontal speed of target, ft/sec 

initial position of target, ft 

present position of target, ft 

present coordinates of missile, ft 


Xjj,y^,z-^ error due to noise in measuring target's position, ft 





V- 


iit/3 

e 


1!he letter 
functions . 


cp with subscripts attached will "be reserved for correlation 
A "bar over any quantity will denote its average value. 


ASSUMPTIONS 
Beam Rider 


We shall assume first of all that the missile it is desired to design 
will be a beam rider. The source of the beam will be considered to be far 
away from the target, so that as the target maneuvers the beam moves only 
parallel to itself, without rotation. 


Target Motion 


It will be assumed here that the target is initially (i.e., at the 
time of firing of the missile) flying with constant speed along some 
straight line. At some time after firing, the target will maneuver. Due 
to the difficulty of maneuvering a large bomber and the short time of 
flight of most missiles, we shall assume that the pilot has time for but 
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one maneuver "before the attack situation is over; for definiteness at this 
time, we take this to he a pitching maneuver 1 in -which the "bomber either 
climbs or dives at some constant acceleration which may or may not be the 
maximum of which it is capable. We assume as an approximation that this 
vertical acceleration leaves the bomber's forward, velocity unimpaired, so 
that the flight path is a parabola rather than a circle. 

Select a coordinate system fixed in space with the X axis parallel 
to the target's initial velocity vector and -the origin at -the initial 
position of the missile. According to what has gone before, we have 


X B » Xq + Vt 


>? 

II 

£ 


( z o > 

23 (z 0 + | a (t-t 0 ) 2 , 

t < t ( 
t > t, 


where t Q denotes the time at which the bomber begins its maneuver. 

The preceding equations represent a whole class of missile -bomber 
combinations. If all the parameters, x©, V, a, etc., are given definite 
values, a particular attack is defined. At this point, we assume that we 
know approximately how far from the missile a target will usually be ini- 
tially, how fast a modem bomber will be going, and how rapidly it can 
accelerate. More precisely, it will be assumed that probability distri- 
butions of the parameters Xq, y Q , z Q , V, and a are known. (As will be 
seen, the entire distributions will not actually be needed, since only 
the mean-square values of these parameters will occur in the work.) 

With Av denoting the average with respect to all parameters, we can 
then compute the following correlation functions. (See ref. 4 for the 
complete definition of these functions.) 

%c B x B ( t > T ) = Av[(x 0 +Vt)(x 0 +VT)] 

= Av(xq 2 ) + Av(xoV)(t+T) + Av(V 2 )tT (l) 

x It will be seen that the argument which follows remains unimpaired 
if the vertical plane In which we have assumed the target maneuvers is 
rotated. Therefore, the analysis which follows will apply equally well 
to situations -where the bomber turns or turns and dives, etc. 
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If the ini tial position of the "bomber is uncorrelated -with its velocity, 
and if, with respect to the missile’s initial position, the bomber is as 
lik ely to be going in one direction as another, we obtain Av(x 0 V) = 0. 
Equation (l) then becomes 

= *o 2 + v2tT 


Similarly, 


<P. 


y B y B 


(t>T) 



To compute cp , it is necessary to assume some probability distri- 

bution for the time t Q at which the target begins its maneuver. 

Although it is possible to make more realistic assumptions, here we shall 
insist that t Q is equally likely to have any value between 0 and T, the 
time at which the missile speed has fallen off so badly that it can no 
longer reach the target. This assumption is not too unrealistic since a 
pilot who is tense will usually maneuver when t Q is small ; the more 
composed will wait longer. 


Knowing these things, we can now write down 
that 0 ^ r 1 t. Then, 


<Pz B z B ’ 


Assume first 


cp z B z B ( t ' T ) = Av[z B (t)z B (T) ] 


i r T 


= -J Av*[zg(t)z B (r) ]dt D 


( 2 ) 


where Av* denotes the average with respect to z Q and a alone. 
Equation (2) then gives 


q W t ' T) 



+ 2 a ("k~"fco) ][ z ° 
+ \ a (t-t 0 ) 2 Jzo| 



= z,. 


■4 


+ 


120T 


(lOt 2 T 3 - 5tT 4 + T 5 ) , 


for 0 ^ t ^ t 
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In deriving this result, we have, 
It follows from equation (2) that 


as "before, assumed that Av(z 0 a) = 0. 
cpz B z B ( t ^) = <Pz B z ;B ( T > t )- Hence, 


U 


w 


9, 


Z B Z B 


(■fc,r) 



a g 

120T 


(lot 3 -^ - 5t 4 T + 1 5 ) , 


for r > t 


We define the correlation function of the homher cpg(t,T ) as the sum 
of the functions cp x ^ x ^,cfy^y^, and Tz^Zg* Thus, 


cp^t,^ = R 2 + V®tr + (I0t 2 r 3 -5tr 4 + r 5 ) 


( 3 ) 


for 0 < t < t. 

It is important to notice that the correlation function ( 3 ) depends 

only on the mean-square values R 2 , , and a 2 . This means that if the 

plane in which the homher 1 s maneuver takes place is not vertical, the 
resulting correlation function will have the same form. Thus, hy choos- 
ing a 2 to he the mean-square homher acceleration over all possible orien- 
tations of the plane in which he maneuvers, we shall he considering the 
more general problem of a homher, flying initially on a straight line 
course, maneuvering in some unknown direction. 


Noise 


The measurements of target position will presumably he made hy radar 
or some other device which will he subject to error. It will he assumed 
that the noise is uncorrelated with the target motion and iB "white.” 
Thus, the autocorrelation of the noise will he given in the x, y, and z 
directions, respectively, hy 


^ X N X N 

T z z 
Z N Z N 


(t,r) 

(t,T) 

(t,T) 


N x 5(t-r) 

NyS(t-T) 


€ 


% 


N z S(t-r) 
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where N x , Ny, and N z are constants, dependent upon the configuration of 
the target, and S(t-r) is the Dirac 5 function (ref. 7)* We define the 
correlation function of the noise as the sum of these three functions: 

cp N (t,r) = MB(t-T) (4) 


■where N = N x + Ny + N z . 


Inputs 


The inputs to the missile consist of the sum of the target motion and 
the noise in each direction. Now, in general, the attack may take place 
from any direction relative to the target, even though certain directions 
may he more probable than others. However, the radar only perceives 
motions as if they occurred in a plane perpendicular to the beam. Thus, 
the information sent the missile is not actual target motion, but the pro- 
jection of this motion onto such a plane. 3h an attack from the beam with 
the missile and the target in the same horizontal plane, these two motions 
coincide. In other situations, the effect of the deviation from a beam 
attack is the apparent reduction of the target velocity and acceleration. 
Thus, we may assume a beam attack by introducing the notion of an apparent 
target whose velocity and acceleration are less than that of the actual 
target. 

Looked at another way, what we are saying is that the procedure of 
optimizing a beam rider may take place as if the plane in which the tar- 
get is moving were orthogonal to the beam. In the final answer, the val- 
ues of V 2 and a 2 which should be used are the apparent values averaged 
over all possible attack directions rather than the true values. 


OPTIMIZATION 

The Minimization Criterion 


The over-all missile-beam system will generally be described by two 
transfer functions, one for the horizontal and one for the vertical direc- 
tion. Although it is theoretically possible to attain a smaller mean- 
square error by actually using two different transfer functions, in this 
report we shall consider only cruciform missiles, so that the transfer 
function is the same in both directions. 

According to the discussion in the preceding section, we may assume 
a beam attack, so that the y axis is parallel to the beam. Let g(t,T) 
denote the impulse response 2 of the missile-beam system. Then, the output 

2 This response depends upon two variables since the system may, in 
general, be time-varying. 
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of the missile will he given by 

->t 


x^Ct) = J g(t,T)[x B (*r) +X N (T)]dT 


5 M("t) = / g(t,T)[z B (T) + Z H (T)]dT 


We desire to minimize the mean-square error 


E 2 = Av[xg(t) - x^(t)] 2 + Av[ z-g(t) - z M (t) ]‘ 


A little algebra yields 

•>t 


E 2 


= cf^Ct^t) - 2 f g(t,r)cp B (t,r)dT + 
o 

J g(t,r) J g(t^ cr)cp B (T,cf)d0 dr + N J g a (t,T)dT , t>0 (5) 


where equation ( 4 ) has been used to give an expression for cpjj. This is 
a minimum (cf. refs. 5 and 8) if and only if 


n t 

<P B ( t > T ) = / g(t,a)qj B (T,0)da + Ng(t,T) 


for 0 < t < t 


(6) 


where cpg is given by equation (3). 

Equation (6) is an integral equation which must be solved for the 
optimum impulse response g(t,r). 


j 


The Optimum 


Although it is a lengthy procedure, the method of reference 4 (or, 
more easily that of ref. 5) can be used to solve equation (6). The opti- 
mum impulse response one obtains is too complicated for there to be much 
advantage in writing it down, especially since a slight approximation 



NACA TN 104-5 


9 


simplifies it enormously. This approximation -will "be discussed in the next 
section. One can say, however, that every thing turns out to be a function 
of the dimensionless time 


t* = At 

where 



( 7 ) 


(8) 


The optimum impulse response g(t,-r) can he used to compute the 
minimum rms error E^y, . Again this function is rather complicated, hut 
figure 1 can he used to find the minimum error as a function of time when 
the optimum missile is fired at a particular target. In order to inter- 
pret figure 1, we note that the mini mum error has been broken up into 
four parts, a part E B which represents error due to initial error at the 
time of firing, a part Ejq- due to noise, a part Ey due to target veloc- 
ity, and a part E B due to target maneuver. The total error is given by 
the square root of the sum of the squares of these four quantities. 

As can be seen from the figure, E B and Ey damp out rather rapidly. 

What remains can be approximated quite accurately by its value at infinity. 
One obtains from the figure that 

e(°°) -- 7% 2 (°°) + e b 2 (°°) 



(9) 


This expression for the minimum error is very important and will be used 
frequently in the sequel. 


AN APPROXIMATION 


Equation (9) can be used as a. check on the performance either of a 
given missile or of one designed by any heuristic process whatever. To 
see whether any missile transfer function, designed by whatever means, 
can be classed as "good" in the sense of this report, one has only to 
check on two points. First he has to see that the rms error, when t = oo, 
is approached rapidly. This is always true of stable systems. Second, he 
may compare the error of his missile with equation (9) to see how far from 
the optimum he may be. 
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Tn -what follows we shall make frequent use of the expressions as ”t 
approaches infinity" and when "t is large." Actually, neither of these 
is quite accurate, since, of course, we shall. never be concerned with 
times greater than T, when the attack is over. What really will be meant 
by large values of t are those values which, while less than T, are 

still far enough from zero that the expression 1 - e"^ may be approximated 
by unity. As can be seen from equation (8), A is not usually very small 
and so "large values of t" may be quite close to zero. 

In this section, we shall utilize the above approximation to reduce 
the formula for the optimum system to manageable proportions. Now, it 
happens that the optimum impulse response depends only on hyperbolic Bines 
and cosines of t* and f* = At. Except near zero time, however, the fol- 
lowing approximations are valid. 


X x 

cosh x = sinh x = — e , 


x > 0 


If these approximations are made in the expression for the optimum impulse 
response and if only fairly large values of t are considered, one obtains 
the simple approximation _ 


g 


;*(t, T ) - + i + | (10) 


where A is given by equation (8) and 


irt/3 1 .43 

e 'I 41 — 


( 11 ) 


if 


It should be noted that since g* is a function of t-T alone, it 
represents a time invariant system. Further, since the real parts of -p 
and +p £ are negative, the system is stable. The transfer function corre- 
sponding to the impulse response (10) is 


2 

+ A 


s + 


2 _ 

A 2 



Writing this in the form 


1+2 3 + 
«n 



(1 + TS) 


(’ 


1+2 



s + 



( 12 ) 


A 


V 
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one obtains the values 


“n 



£ 


n 



u d = A 



1 

2 


T 


1 

A 


■where A is given by equation (8). 

The error corresponding to equation (10) can be computed by substi- 
tuting (10) into (5). The exact procedure for doing this for any time 
invariant system will be indicated later ohj it suffices here to say that 
the error corresponding to g(t,-r) for large t is given again by for- 
mula (9). Thus , although the error of this approximate system may be 
larger than the minimum initially, this transient phase soon dies out and 
the error for large t is compromised not at all. 


ACCELERATION DEMANDED OF THE MISSHE 


We have thus far specified an impulse response (10), corresponding 
to the transfer function (12), which has an associated error (9). However, 
we have no guarantee that in order to achieve the performance indicated 
by (9) the missile may not be called upon to attempt impossible exerti ons . 
The rms acceleration of the missile, for example, may be required to be 
so great that the controls of the actual missile will be at the stops an 
of the time, thus ma k ing the linear analysis used here invalid. The sys- 
tem may, of course, be limited in other ways, but we shall consider only 
acceleration limiting as typical of the sort of ' situation which arises, 
and, in addition, as the most important type of limiting. 

Now, whatever the inputs and whatever the impulse response g(t,-r) 
may be, the mean-square acceleration demanded of the missile is 

®M 2 = Av “£2 f s(b, T ) [x-g( T ) + Xjj(t) ]dTj>- + 


s(b,T)[z B (T) +z K (T)]dt| 


(13) 
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Define 


Xj(t) = Xg(t) + Xjj(t) 
Z I (t) = ^(t) + Zj^t) 


(3A) 


Substituting these definitions into (13) and differentiating, one obtains 




, 2 * Av 


■|g(t,t)x x (t) + [2gj. (t, t) + g T (t,t) JxjCt) + J SttCt^TjxjCrJdrj- 

Av |g(t,t)z I (t) + teg^tjt) + g T (t,t) IzjCt) + J ^(t^TjzjCrJdT^ 


■where the subscripts on g denote partial differentiation. 

Now let cp- . ( t . t) denote the autocorrelation function of x T , 
T Xj*X*j“ 

(p».^ y (t,r) the cross -correlation function of Xj with Xj, etc. Then, 
squaring and. averaging, one obtains 


j^2g t (t,t)+6 T (t,t)j 2 J^q^ i (t,t)+<5) ! , iZi (t,t)J + 2^ g(t, T )|g(t,t)J^q) jtiXi (t,T)+<pj iI . l (t,T)J + 

[2g t (t,t)+g T (t,t)J<p x ^(t,T)^ ZiZi (t,T)jj-aT+^ g fct (t, T )|' Sb t (t,o)^<p XiXi (T,a)+<p ZiZi (T,o)jda Or 

(15) 


Consider now the first bracketed, term in this equation. It contains, 
according to equations (14), a term involving the noise. This term is 
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+<p w t ' t) 


AvtXjjCtJXjjCT) + Z H (t)Zjj(T) ] 


T=t 


a 8 

St St 


Av[x N (t)x H -( T ) + z N (t)z N ( T ) ] 


T=t 


S 2 

St St 


%(^t) 


T=t 


(16) 


Now, if the noise is white, cp^ is given "by equation ( 4 ) and so (l 6 ), 

which is to "be evaluated at t = t, is infinite. Even if the strictly 
unrealizable assumption that the noise is white is not made, it can be 
expected that this term will be unduly large. In order to e limina te it 
then, we set, according to equation ( 15 ), 

g(t,t) = 0 (17) 


Similarly, the third bracketed term in ( 15 ) gives 

2 g t (t,t) + &j-(t,t) = 0 ( 18 ) 

We now specify (17) and (18) as inviolable conditions which any satisfac- 
tory missile must satisfy in order that the missile acceleration which is 
called for be not too large. In this case, ( 15 ) gives for the acceleration 
of the missile 



“/ «bt ^t^^x-^*) + < Pz I z I (^ cr ) 


dcr dT 


(19) 


Let us see if the missile defined by the impulse response g* of 
equation (10) satisfies conditions ( 17 ) and (l 8 ). We have from (10 ) tha t 


g*(t,t) = ^ 

= j (5 + t - P 2 ) 


= 2A 
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since' \i = 3 . The constant A 0, meaning that condition (17) is not 

satisfied. - Thus, we may conclude that the acceleration demanded of the 
system described by ( 10 ) will be so large that the system is unsatisfactory. 


OPTIMIZATION IN THE PRESENCE OF CONDITIONS (17) AND (l8) 


We have seen that the system which has been derived turns out to be 
unsatisfactory from the point of view of acceleration limiting. It appears 
that we might proceed using the ideas of reference 9 , minimizing not -tile 
error E 2 , but the quantity 

E 2 + f (20) 


where 7 is a constant. The idea here is that when 7 is large, effec- 
tively we would be minimizing a^ 2 while, when 7 is small, E 2 would be 

minimized. Consequently, the thought is that there might be a value of 7 
which while keeping E 2 near the minimum value (9) still does not allow 

a M 2 ' to get too large. However, in our case, the quantity (20) has no 

minimum, as can be shown by a straightforward (though lengthy) computation. 
Hence, we must proceed otherwise. 

What we shall do is utilize a hint given us by the form of g(t,-r) 
and actually compute the error using formula (5) . It will turn out that 
conditions necessary In order that E 2 be small will then become evident 
by inspection. 

The hint is that g(t,r) represents a time -invariant system. Thus, 
if we decide to look only at time -invariant systems, we may expect, in 
view of the fact that (12) is time -invariant, to be able to come upon a 
fairly good result. 


Optimization 


Now, the most general time -Invariant system has an impulse response 
which is a sum of exponentials, provided only that the roots of the 
characteristic equation are distinct. 3 Hence, we set 


s(t,T) = g(t-r) = E Bn e"* n ^ t “ T ^ 


( 21 ) 


^ven if they are not, the form of the impulse response is a limit 
of a sum of exponentials as the roots become coincident. 
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We may now substitute (21) into the expression ( 5 ) for the error. 
This procedure is simplified by breaking the error up into four parts, 
as in figure 1. Thus, we define, according to the expressions ( 3 ) 
and ( 4 ) for the correlation functions cpg and cp^-, 


E R 2 (t) S R 2 


£ 1 - 2 ^ g(t,T)dr +£ g(t,T)J^ g(t,cr)dcr ItJ 

= i 2 [l - / tg(t?T)dT ] 


Ey 2 (t) = 


V 2 j^t 2 - 2t/ rg( t, t ) dT +£ g(t,T)/ g(t,a)TCT da ItJ 

= V 2 jjt - J Tg(t,T)drJ 


>( 22 ) 


% 2 (t) 


N f g 2 (t,r)dr 
o 


E/(t) s -g- j^6tF - 2j rt g(t,r)(l0t 2 T 3 -5tT 4 + T 5 )dT 



g(t, a) (10t 2 c 3 - 5tct 4 + a 5 )dcr 
g(t,cr) (lOtV - 5 t 4 ct + t 5 ) dcr 


dT + 


In equation (21), g(t,T) is a function g(t-T) of the difference of t 
and. t. This means that formulas (22) can be written 



1 6 


NACA TN 4145 


*fc "1^ 

E R 2 (t)=p[l -J e(T)drJ 


-b 1 'l 2 

s/(t) -/ eW iT J + ^ Tg(T)dT| 


%■ 


pt 

: (t) = Hf / ^(t)^ 


2 r t g(T)^ T g(«H ( t- T ) s - 5(t-T) 4 (t-o) + ia(t-T) 3 (t- ff ) 2 ]a.a dr} 


ph 

I n (t) = J T^girl&T , n = °A> • • ’ 

o 

Then. 

% 2 (t) = R 5 [1 - 

■\ 2 

Ey 2 (t) = V 2 «jt[l - IoC^n + I i( t )j 

pt 

E w 2 (t) = N J f?{r)lT 
o 

%«(*> -5^{ 6t “ti-ioWl 2 ^ti-io(*)3^W ' 

^[A-I^t)} I 2 (t)-2I 1 2 («)J-6ot 2 I 1 ( 1 >) I 2 ( t ) + 3 0tI 2 2 { t ) + 

£ r t g (T ) Wl-Io ( T)}+^ 4:l; 1 . (T) - :L0T3l 2 (T) ] aT } 

J o 
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Now, if we are to set g(t) equal, as in equation (21), to a sum of 
exponentials, and if the system is to he stable, all of the quantities 
I n (t) will approach definite limits as t Consider, then, the t 5 

term in the expression for E B 2 . Since Io(t) approaches a limi t, this 
whole term will approach infinity as t->- co unless l-I 0 (eo) = 0. This 
gives us one relation which g(t) must satisfy: 

I 0 (“) = 1 (24) 

Similarly, the t 3 term in Eg 2 gives 


while the t term gives 


= 0 


I 2 («) = 0 


When these conditions acre satisfied, we can say that 

% 2 (°°) =0 \ 

E v 2 (<x>) =0 I 

En^oo) g 2 (r)dT ^ ( 2 

o 

E b 2 (oo) = ^^J' s(t)/t 5 [1-I 0 (t)3 + 5t 4 I 1 (t) - 10t 3 I 2 (t) 1 dx 

o'- J 

The acceleration (19) of the missile can also he broken into four 
parts. Recalling that if g(t,T) = g(t-T), then g t . + .(t,T) = g(t-T), we 


see that 


n t pt 

a R 2 = R 2 / g(t-T)/ g(t-cr)dcr dT 

o o 

= R 2 j ' g(T)drJ 
= R 2 J^fe(t) - g(0)J 


= R 2 g 2 (t) 
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by condition (l8), which 

reads, in this time -invariant case. 



0 

11 

o' 

•5) 

(29) 

Similarly, (17) becomes 

g(0) = 0 

(30) 

and so one can write 

a v 2 = V 2 g 2 (t) 

(31) 

by using condition (30) . 

Also, 



a N 2 = N ^ g 2 (T)dT 

(32) 


o 


Now, set 

J n ( t ) = / T n g(r)dr 

J o 

= t n g(t) -nt n " l g(t) +n(n-l)l n . 2 (t) (33) 


Then, just as before, 




g(n)[ (t-T) 5 -5(t-T) 4 (t-ff)+10(t-T) 3 (t-a) 2 ]da dr 


= -§£. X 

6ot L 


6t 5 J o 2 (t)-30t 4 J o (t)J 1 (t)+20t 3 [J o (t)J 2 (t)+2J :L 2 (t) ] - 

t 


6ot 2 J 1 (t)J 2 (t)+30tJ 2 2 (t)-2^' g ( t) [ 6t 3 g ( t) -15t 4 g ( t ) +20t 3 Io ( t ) ]drj- 


From (33), however, we see that If the system is stable (so that 

g(oo) = g(oo) =0), then J 0 (°°) = J-^oo) = 0. Also, using conditions (2k), 

(25), and (26), we obtain J 2 (co) = 2. Therefore, as ~, 


a B 2 (t) 


6ot 




g(f) [6T 5 g(T)-15T 4 g(T)+20T 3 I o (T) ]dT 


} 


m 
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Now, let g(t) be given by (21). Conditions (24), (25), (26), (29), 
and ( 30 ) then become 



Notice that the first three of equations ( 35 ) are conditions for the error 
not to be too large, while the last two are conditions for the acceleration 
to be finite. Finally, we have 

E 2 (oo) = Eft 2 (00) + Eg 2 (») 




where 4 A is given by ( 8 ) and the sum runs over al 1 m and n. Since we 
still are interested in reducing E 2 as far as possible, the problem now 
becomes one of minimizing ( 36 ) subject to the conditions ( 35 ) . 

Equation (21) can be substituted into ( 28 ), (31), (32), and ( 34 ) to 
yield an expression for the acceleration of the missile. Assuming sta- 
bility, we see that and a v 2 quickly approach zero. On the other 

hand, 

a^- 2 (t) ~ ajq- 2 + a-g 2 


= a 2 i + N 


vV 2 + ^ 


for large t. ~ - ■ 


4 An alternative formulation of formulas ( 36 ) and ( 37 ) is given in 
the appendix. 
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It is very important to notice that there are infinitely many systems 
satisfying conditions (35) . In fact, even if a characteristic equation 
(of sufficiently high order, see below) is given for the determination of 
the roots A n , there still exists a system (perhaps many such) satisfy- 
ing (35). To see this, consider the 7^'s as fixed. Equations (35) then 
give relations among the a^'s. These, however, determine only the numer- 
ator of the transfer function. Hence, we can say that given a denominator 
one can find a corresponding numerator such that the system satisfies (35)* 

This result shows that not only can one use the ideas discussed here 
to design a missile but that one can very easily add a compensating net- 
work to an existing system so that the new over-all system will satisfy 
the conditions (35) . This follows from the fact that equations (35) only 
give relations among the a^'s. Further, since the Vs are still arbi- 
trary, they remain as parameters with respect to which the error ( 3 6 ) can 
be minimized. 


The Transfer Function 


Conditions (35) may easily be interpreted as conditions on the 
transfer function G(s) corresponding to the impulse response g(t) . In 
fact, since 

g(t) = Za^e ^ 


we have 



Sn 

s+A n 


Consequently, the first three conditions (35) become 

G(0) = 1 ' 

G'(0) = 0 
G"(0) = 0 , 


(38) 


Consider next the last two of equations (35) • They say that 
S 1 (0) = g(0) = 0. As is well known, however, this can only be so if the 
order of the denominator of g exceeds the order of the numerator by at 
least 3. Thus, we need only add the relation 

(order of denominator of G) - (order of numerator of G) ^ 3 (39) 

to equations ( 38 ) to assure ourselves that all of equations (35) 
satisfied. 


are 
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Now, suppose 


G(s) 


P(s) 

Q(s) 


■where P and Q are polynomials of order m and n, respectively. Equa- 
tion ( 39 ) then reads 

n - m > 3 (40) 

■while equations ( 38 ) "become 

P(0) = Q(0) 

P*(0) = Q'(0) 

P"(0) = Q" (0) 


Thus, equations ( 38 ) are satisfied if the first three terms (i.e., the 
constant, linear and quadratic terms) in P(s) and Q(s) are equal. We 
may, by dividing through by the constant, assume it to be unity and say 
that conditions (35) will be satisfied if G(s) has the form 


G(s) 


m-3 

l-K^s+OaS^s 3 / p v s v 

l+a 1 s+a a s 2 +s 3 \ 
v=o 


(*H) 


and (40) holds. 


Note that it follows from (40) and (4l) that the order of the system 
is at least 5> since obviously, in view of (40), the lowest order system 
will occur if m = 3 and p Q = 0 - in which case, again by (40), n must 
equal 5 * a t least. 


SOME OPTIMUM TRANSFER FUNCTIONS 


Tn this section, we wish actually to display some systems satisfying 
conditions ( 35 ) and minimizing ( 36 ), in order to show the effect on the 
missile acceleration. 

As shown in the prece din g section, any such system must be of the 
fifth order, at least. For simplicity, we shall only consider systems 
with this minimum order. Thus, according to equation (4l) and the remarks 
following, we write 
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1+a. s+a^s 2 

8(b) 

l+a 1 s+a s s 2 +a 3 s 3 +a 4 s 4 +05 s 5 


(42) 


We know that the impulse response corresponding to (42) satisfies (35) If 
a 5 ^ 0. Notice that, except for the further condition of stability, the 
five constants a are free and can be used to minimize (36) . It is pos- 
sible that with so much freedom, the error (36) - even when the constants 
& n are subject to conditions (35) - can be reduced to the absolute minimum 
error (9). However, the minimization of (36) in this way is extremely com- 
plicated and the problem probably can only be solved by use of a method 
such as steepest descent (ref. 10) on a digital computer. Here, we shall 
solve a much less ambitious problem; however, we shall show that in order 
to reduce the missile acceleration very greatly, a large penalty in 
increased error need not be paid. 


What we shall actually consider are special cases of (42) having the 

form 



(43) 


where, as before, A is given by (8). The transfer functions (43) depend 
on a parameter p. Notice that as p-> 00, (43) reduces to the optimum 
transfer function (12). Thus, for large P at least, one may expect 
the error corresponding to (43) to be somewhere in the neighborhood of 
the minimum error (9) • 

The act ual error of the system (43) can be computed with the aid 
of (36). The computation, though long, is not impossible, and the final 
result is given in normalized form in figure 2. Note that as predicted, 
the curve is fairly flat, except at the very left-hand end. 

In order to understand what significance the magnitude of p has, 
it is necessary to look at the mean-square acceleration of the missile. 
This was computed, using (37) with t = T (the worst case). The result is 
shown in normalized form in figure 3 (the source of the parameter used as 
ordinate can be seen in the appendix) . Figure 3 also contains the curve 
of figure 2. ' r - '• r ' i - i; ; - - 
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Figure 3 represents the climax of. all our efforts. In order to use 
the figure, one would have to decide what the rms acceleration capabilities 
of the missile were to he. Using this knowledge, one finds the appropriate 

( a M 2 “ a2 )/NA 5 and reads off the corresponding value of abscissa p. This 
value of p is then used in equation (43) to give the optimum transfer 
function. The error corresponding to this system is then read at the same 
value of p. An example of this is given in a later section. 


The Error 


The error of the system (li-3 ) is plotted against the dimensionless 
missile acceleration parameter A = (a^ 2 - a 2 )/NA 5 in figure 4. Notice 

that the error E is dete rmin ed by only two things, the minimum error 
E m j_ n and this parameter A. An advantage of the present applications 
of the method of "optimization" is that the parameters of which E is 
a function are determined analytically. 

Note further that over what might be considered some reasonable range 

of A - say, those values corresponding to the 13.nd.ts 4 < 7, 

a condition which with a "reasonable" choice of a 2 , N, and T is the range 
bounded by the dotted lines 3n figure 4 - the slope changes slowly. One 
of the principal advantages of the present outlook is that the curve of 
figure 4 gives the error over the practical range of A uniformly. 


AN EXAMPLE 


As an example of the use of figure 3> consider a target having the 
following characteristics 5 



N = 100 ft 2 sec 


5 The value of N given here is a reasonable one and is consistent 
with that used in earlier reports. In reference 1, for instance, the 
value of "noise magnitude" which was used was 15 ft 2 / radian/sec ; since 
our definition of N differs from that used in reference 1 by a factor 
of 2 3t we have set N = 2«(15) ~ 100. 
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Let us assume that the time of flight T is 10 seconds, and that the rms 
value of missile acceleration will he 7s' s > 80 that 


= 225 ft/sec 2 


We then have 



= 1.00 


Also, 


2 - a£ 
NA 5 


498 


Hence, reading along the dotted lines in figure 3> we see that we must 
choose P = 8.6, which gives a value of the error of 

E = 1.19 


and, since 


E min = ^2HA 


= 14.1 ft 


we have 

E = 16.8 ft (44) 

With p = 8.6 and A = 1.00, the transfer function (43) becomes 


G(s) 


1+2.23b+2.49b 2 

(1 +b) ( l+s+s 2 ) (1+0. 23 b+ 0. 027s 2 ) 


m 


In the attack situation described, the system having the transfer function 
(45) will have the error (44). 



HA.CA TN 4114-5 


25 


CONCLUDING REMARKS 


It has heen shown how a team rider missile system can "be designed in 
an optimum way even when a target maneuver is a strictly nonstationary one. 
The maneuver was chosen here to be a step in target acceleration, the step 
being assumed to occur with equal likelihood anywhere in a finite interval 
of time. A simple formula, was derived for the minimum error in these cir- 
cumstances; this formula can be used as a criterion of merit against which 
the value of any beam rider, designed by whatever means, can be measured. 

The optimum control system was shown to be unsatisfactory Bince the . 
accelerations demanded of the missile were unduly large. Therefore, it 
was necessary to design it in a different way, taking into account the 
limited acceleration available to any real missile, and using the minimum 
error formula to decide how good the new system was. Two things resulted 
from this analysis. The first was the conclusion that the transfer func- 
tion of any missile guidance system must be of a certain form (given by 
eqs. (40) and (4l)) if either the missile acceleration or the miss distance 
is not to increase beyond all bounds. 

The second result was presented in figure 3 and equation (43). Fig- 
ure 3 may be used as a design chart to determine a satisfactory system 
transfer function in any given situation. Thus, having decided just how 
much acceleration the desired missile would be able to withstand, the 
designer can use figure 3 to determine a corresponding value of a cer- 
tain parameter (3. The transfer function is then determined by means of 
equation (43). The error of a missile system with this transfer function 
may then be read off figure 3 • 

Finally, the curve of error versus a dimensionless acceleration 
parameter - which is the only thing which determines the error - is 
given. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
Sept. 30, 1957 
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APPENDIX 

ALTERNATIVE FORMULATION OF EQUATIONS (36) AND (37) 


In order to evaluate any system, it is necessary to calculate the 
error corresponding to it (of course) and its rms acceleration. With 
the proviso that conditions (29) and (30) hold, formulas (36) and (37) 
of the text can he used for this purpose. However, it is simpler compu- 
tationally to consider a transformed version of these formulas . This 
version is described 1 below. 


E‘ 


; (oo) =■ N 


m.n 


N 


(36). 

We have 


a m a n ( 

(l + * \ 


Affl+A n | 



a m a n / 
0 


-Ota 

e 

r . 

v 0 

/ 


V_fn 

k ZjA n 3 


-An) t 


dt 


dt 


and in view of equation (21), 

00 

= N J' [g 2 (t) + A®h 2 (t) ]dt 


where 



Now, by a well-known formula for Laplace transforms (ref. 11 ) 
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where G(s) is the Laplace transform, of g(t) . Also, again "by (21), 

rpo 

kfa) = (t-r) 2 g(T)dT 


Consequently, the transform of h(t) is 


H(s, -/*«(*) [£-£♦£*«£] 


dt 


= G"(0) G'(0) GM _ G(s) 

2s B 2 s 3 s 3 


However, as shown on page 20, conditions ( 35 ) entail (38), so that if 
conditions (35) are valid for the system under consideration. 


H(s) 


l-G(s) 


Also, 


/ 


, 10 ° 


h 2 (t)dt = J H(s)H(-s)ds 

-ioo 


Therefore, finally. 


E 2 (oo) = JL 
2jcL 


[tow-.)-* **;' 


-ioo 


N 

2iti 


— — ico 


-ioo 


G(s)o(- B )as 


-ioo 


l-G(s) l-G(-s) 
s 3 s 3 


ds (Al) 


hy virtue of the definition (8) of A. 
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In a similar way, it can "be shown that equation (37) is equivalent 
to 

^?(*) - ^ J 1 ” U(s) B (-s, - x 6 hghil j as (A2) 

-ico L 




l-G(s) l-G(-s) 
s s 




provided again that ( 29 ), ( 30 ), and the first of equations ( 38 ) hold. 


The in tegrals occurring in equations (Al) and (A3) can he evaluated 
immediately hy means of the tables E.2-1 in Appendix E of reference 12. 


One further point which might he made is this. It is easy to show 
hy use of equations (A2) where the parameter (a^-a 2 )/]?! 5 came from. 
Notice that in equation (43) that G(s) is really a function of s/X, so 
that G(Xs) does not involve X explicitly at all. Hence, writing s =«Xp 

In (A2) and abbreviating a^(T) to a^s, we obtain 





■loo 


p 4 g(Xp)g(-xp) 


l-G(Xp) l-G(-Xp) ' 
P P 


dp 


or 



1 

2 id 



p 4 G(Xp)G(-Xp) 


l-G(Xp) l-G(-Xp)"L_ 

dp 

P P 3 


where the right-hand side is in a normalized form and, contrary to 
appearances, does not contain any of the parameters a 5 ”, N, T, or X. 
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